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Report 

We  achieved  considerable  success  in  developing  new  kinds  of  density 
functionals  and  new  density  functionals  for  catalysis.  We  also  made  several 
successful  applications  of  density  functional  theory. 

As  a  long-standing  puzzle,  experimental  observations  reveal  faster 
organophosphine  dissociation  in  the  olefin  metathesis  by  Grubbs’  s  first- 
generation  precatalyst  (Gen  I)  than  by  the  second-generation  precatalyst  (Gen  II), 
but  Gen  I  shows  less  catalytic  activity.  We  showed  by  electronic  structure 
calculations  with  the  M06-L  density  functional  that  carbene  rotamer  energetic 
effects  are  responsible  for  the  inverse  relation  between  organophosphine 
dissociation  rate  and  catalytic  activity.  The  carbene  rotamer  acts  as  a  toggle 
switch,  triggering  the  dissociative  mechanism  that  produces  the  active  catalyst. 
The  slower  catalyst  production  in  Gen  II  as  compared  to  Gen  I  is  not  a  pure 
electronic  effect  but  results  from  rotameric  coupling  to  the  dissociation  coordinate 
speeding  up  Gen  I  dissociation  more  than  Gen  II  dissociation.  If 
organophosphine  dissociation  were  to  occur  with  fixed  rotamer  orientation,  Gen  II 
would  be  produced  faster  than  Gen  I,  as  originally  expected.  The  rotameric 
energetics  also  contributes  to  the  higher  catalytic  activity  of  the  Gen  II  catalyst. 

For  molecules  containing  the  fourth-period  element  arsenic,  we  test  (i,  ii)  the 
accuracy  of  all-electron  (AE)  basis  sets  from  the  def2-x  ZVP  and  ma-x  ZVP 
series  (where  x  Z  is  S,  TZ,  or  QZ),  (iii)  the  accuracy  of  the  6-31 1 G  series  of  AE 
basis  sets  with  additional  polarization  and  diff  use  functions,  and  (iv)  the 
performance  of  effective  core  potentials  (ECPs).  The  first  set  of  tests  involves 
basis-set  convergence  studies  with  eleven  density  functionals  for  five  cases: 
equilibrium  dissociation  energy  (De  )  of  As2  ,  vertical  ionization  potential  (VIP)  of 
As2  ,  IP  of  As,  acid  dissociation  of  FI3  As04  ,  and  De  of  FeAs.  A  second  set  of 
tests  involves  the  same  kinds  of  basis-set  convergence  studies  for  the  VIP  and 
De  values  of  As3  and  As4  clusters.  Both  relativistic  and  nonrelativistic 
calculations  are  considered,  including  in  each  case  both  AE  calculations  and 
calculations  with  ECPs.  Convergence  and  accuracy  are  assessed  by  comparing 
to  relativistic  AE  calculations  with  the  cc-pV5Z-DK  or  ma-cc-pV5Z-DK  basis  and 
to  nonrelativistic  AE  calculations  with  the  cc-pV5Z  or  ma-cc-pV5Z  basis.  The 
primary  objective  of  this  study  is  to  evaluate  the  abilities  of  ECPs  with  both  their 
recommended  basis  sets  and  other  basis  sets  to  reproduce  the  results  of  all¬ 
electron  relativistic  calculations.  The  performance  of  the  def2  and  ma  series 
basis  sets  is  consistent  with  their  sizes,  and  quadruple-^  basis  sets  are  the  best. 


The  def2-TZVP  basis  set  performs  better  than  most  of  the  6-311 G  series  basis 
sets,  which  are  the  most  commonly  used  basis  sets  in  the  previous  studies  of 
arsenic  compounds.  However,  relativistic  def2-TZVP  calculations  are  not 
recommended.  The  large-core  ECPs,  which  are  the  only  available  ECPs  for 
arsenic  in  the  popular  Gaussian  program,  have  average  errors  of  9_  12  kcal/mol 
for  the  arsenic  systems  studied;  therefore,  these  ECPs  are  not  recommended. 
The  triple-^  small-core  relativistic  ECP  (RECP)  basis  set  cc-pVTZ-PP  is  found  to 
have  performance  better  than  that  of  the  def2-TZVP  basis  set,  and  it  is  highly 
recommended  for  arsenic-containing  systems.  The  double-^  RECP  basis  set  ma- 
sc-SVP  is  recommended  for  large  arsenic  systems.  The  double-^  RECP  basis 
set  ma-sc-SVP  is  recommended  for  large  arsenic  systems  for  which  the  def2- 
TZVP  and  cc-pVTZ-PP  basis  sets  are  unaffordable,  if  a  basis-set  error  of  ~  2 
kcal/mol  can  be  tolerated. 

The  Gibbs  energy  of  reaction  of  oxidative  addition  of  ammonia  to  an  iridium 
complex  in  diethyl  ether  was  calculated  by  seven  density  functional  methods,  in 
particular  B3LYP,  PBE,  CAMB3LYP,  M05,  M06,  M06-L,  and  oB97X.  The 
calculated  free  energies,  based  on  geometry  optimization  and  frequency 
calculations  in  both  the  gas  phase  and  solution,  were  compared  with  the 
experimental  result,  -1 .3  kcal/mol,  obtained  by  Hartwig  and  coworkers.  The  M06- 
L  method  gives  the  best  result:  -1 .4  kcal/mol. 

We  employed  noncollinear  density  functional  theory  to  show  that  the  low-spin 
state  of  Mn3  in  a  model  of  the  oxygen-evolving  complex  of  photosystem  II  avoids 
frustrated  spin  coupling  by  adopting  a  noncollinear  arrangement  of  spins,  thereby 
lowering  the  energy  by  7  kcal/mol.  The  high-spin  state  also  has  noncollinear 
spins.  The  optimum  self-consistent  field  solutions  for  this  multinuclear 
oxomanganese  complex  correspond  to  states  that  cannot  be  described  by  the 
unrestricted  Slater  determinants  used  in  Kohn-Sham  collinear  density  functional 
methods.  This  kind  of  spin  coupling  can  be  important  in  many  open-shell 
systems,  and  the  conventional  collinear  spin  interpretation  of  chemical  bonding  in 
such  systems  should  be  viewed  with  caution. 

The  Minnesota  family  of  exchange-correlation  functionals,  which  consists  of  meta 
generalized  gradient  approximations  (meta-GGAs)  and  global-hybrid  meta- 
GGAs,  has  been  successful  for  density  functional  calculations  of  molecular 
structure,  properties,  and  thermochemistry,  kinetics,  noncovalent  interactions, 
and  spectroscopy.  Here,  we  generalize  the  functional  form  by  using  range- 
separated  hybrid  meta-GGA  exchange.  We  optimized  a  functional,  called  Mil, 
with  the  new  form  against  a  broad  database  of  energetic  chemical  properties  and 
compare  its  performance  to  that  of  several  other  functionals,  including  previous 
Minnesota  functionals.  We  required  the  percentage  of  Hartree-Fock  exchange  to 
be  100  at  large  interelectronic  distance,  and  we  found  an  optimum  percentage  of 
42.8  at  short  range.  Mil  has  good  across-the-board  performance  and  the 
smallest  mean  unsigned  error  over  the  whole  test  set  of  332  data;  it  has 
especially  good  performance  for  main-group  atomization  energies,  proton 


affinities,  electron  affinities,  alkyl  bond  dissociation  energies,  barrier  heights, 
noncovalent  interaction  energies,  and  charge-transfer  electronic  excitation. 

The  performance  of  popular  Hartree-Fock-based  effective  core  potentials  in 
Hartree-Fock  and  density  functional  calculations  of  3d  transition  metals  was 
evaluated  by  basis-set  convergence  studies  for  ten  cases:  the  equilibrium  bond 
dissociation  energy  (De)  for  dissociation  of  ground-state  Ti2  to  ground  and  excited 
atoms,  the  ground-state  dissociation  energies  of  FeO,  Cu2,  ScFI,  TiFI,  Sc2,  Fe2, 
and  TiV"^,  and  the  first  excitation  energy  of  Ti  atom.  Each  case  was  studied  with 
11  or  13  density  functionals.  For  comparison,  the  accuracy  of  the  all-electron 
def2-TZVP  basis  set  is  tested  with  both  relativistic  and  nonrelativistic  treatments. 
Convergence  and  accuracy  are  assessed  by  comparing  to  relativistic  all-electron 
calculations  with  a  nearly  complete  relativistic  basis  set  (NCBS-DK,  which 
denotes  the  cc-pV5Z-DK  basis  set  for  3d  metals  and  hydrogen  and  the  ma-cc- 
pV5Z-DK  basis  set  for  oxygen)  and  to  nonrelativistic  all-electron  calculations  with 
a  nearly  complete  nonrelativistic  basis  set  (NCBS-NR,  which  denotes  the  cc- 
pV5Z  basis  set  for  3d  metals  and  hydrogen  and  the  ma-cc-pV5Z  basis  set  for 
oxygen).  As  compared  to  NCBS-DK  results,  all  ECP  calculations  perform  worse 
than  def2-TZVP  all-electron  relativistic  calculations  when  averaged  over  all  130 
data  (13  functionals  and  ten  test  cases).  The  compact  effective  potential  (CEP) 
relativistic  effective  core  potential  (RECP)  combined  with  a  valence  basis  set 
developed  for  the  many-electron  Dirac_  Fock  (MDF10)  RECP  performs  best  in 
effective  core  potential  calculations  and  has  an  average  basis-set 
incompleteness  error  of  3.7  kcal/mol,  which  is  much  larger  than  that  (0.9 
kcal/mol)  of  def2-  TZVP  relativistic  all-electron  results.  Hence,  the  def2-TZVP 
relativistic  all-electron  calculations  are  recommended  for  accurate  DFT 
calculations  on  3d  transition  metals.  In  addition  to  our  general  findings,  we 
observed  that  all  kinds  of  density  functionals  do  not  show  the  same  trends.  For 
example,  when  ECPs  are  used  with  hybrid  functionals,  which  sometimes  are  not 
recommended  for  calculations  of  transition  metal  systems,  they  are  found  to 
perform  better  at  achieving  the  basis-set  limit  than  when  used  with  local 
functionals  and  meta-GGA  functionals.  The  most  successful  combination  of 
RECP  and  basis  set  has  a  basis-set  incompleteness  error  of  1 .7-2.4  kcal/mol  for 
hybrid  generalized  gradient  approximations,  which  is  smaller  than  that  of 
nonrelativistic  NCBS  calculations  (whose  average  basis-set  incompleteness  error 
for  hybrid  functionals  is  2.7-2. 9  kcal/mol).  The  average  basis-set  incompleteness. 

We  reported  a  new  local  exchange-correlation  energy  functional  that  has 
significantly  improved  across-the-board  performance,  including  main-group  and 
transition  metal  chemistry  and  solid-state  physics,  especially  atomization 
energies,  ionization  potentials,  barrier  heights,  noncovalent  interactions, 
isomerization  energies  of  large  molecules,  and  solid-state  lattice  constants  and 
cohesive  energies. 


The  ability  of  density  functional  theory  to  predict  the  relative  energies  of  different 
spin  states,  especially  for  systems  containing  transition  metal  atoms,  is  of  great 


importance  for  many  applications.  In  order  to  sort  out  the  key  factors  determining 
accuracy,  we  compared  the  predictions  of  60  density  functional  approximations 
of  10  different  types  [local  spin  density  approximation,  generalized  gradient 
approximation  (GGA),  nonseparable  gradient  approximation  (NGA),  global-hybrid 
GGA,  range-separated  hybrid  GGA,  range-separated  hybrid  GGA  plus  molecular 
mechanics,  meta-GGA,  meta-NGA,  global  hybrid  meta-GGA,  and  range- 
separated  hybrid  meta-GGA]  for  their  ability  to  represent  the  spin-flip  transitions 
of  all  4d  transition  metal  atoms  of  groups  3-  10  (Y  through  Pd)  and  their  singly 
positive  cations.  We  considered  all  16  excitation  energies  connecting  the  ground 
states  (of  the  neutral  atoms  and  the  cations)  to  their  first  excited  states  of 
different  multiplicities,  and  we  also  consider  all  eight  ionization  potentials.  We 
also  tested  the  Hartree-Fock  method.  All  density  functional  and  Hartree-Fock 
calculations  are  converged  to  a  stable  solution,  in  which  the  spatial  symmetry  is 
allowed  to  be  completely  broken  to  achieve  the  lowest  possible  energy  solution. 
By  analyzing  the  fractional  subshell  occupancies  and  spin  contaminations,  we 
were  able  to  sort  out  the  effects  of  s  orbital  vs.  d  orbital  bias  and  high-spin  vs. 
low-spin  bias.  A  reliable  functional  should  have  little  or  no  bias  of  either  type 
rather  than  succeeding  for  a  limited  subset  of  cases  by  cancellation  of  errors.  We 
found  that  the  widely  used  correlations  of  spin  splittings  to  percentage  of 
Hartree-  Fock  exchange  are  not  borne  out  by  the  data,  and  the  correlation 
functionals  also  play  a  significant  role.  We  eventually  concluded  that  SOGGA1 1- 
X,  B1LYP,  B3V5LYP,  and  MPW3LYP  are  the  most  consistently  reliable 
functionals  for  balanced  treatments  of  4d  transition  metal  atoms  and  their 
cations. 

A  notorious  failing  of  approximate  exchange-correlation  functionals  when  applied 
to  problems  involving  catalysis  has  been  the  inability  of  most  local  functionals  to 
predict  the  correct  adsorption  site  for  CO  on  metal  surfaces  or  to  simultaneously 
predict  accurate  surface  formation  energies  and  adsorption  energies  for 
transition  metals.  By  adding  the  kinetic  energy  density  t  to  the  density  functional, 
the  revTPSS  density  functional  was  shown  recently  to  achieve  a  balanced 
description  of  surface  energies  and  adsorption  energies.  We  showed  that  the 
older  M06-L  density  functional,  also  containing  t,  provides  improved  surface 
formation  energies  and  CO  adsorption  energies  over  revTPSS  for  five  transition 
metals  and  correctly  predicts  the  on-top/hollow  site  adsorption  preferences  for 
four  of  the  five  metals,  which  was  not  achieved  by  most  other  local  functionals. 
Because  M06-L  was  entirely  designed  on  the  basis  of  atomic  and  molecular 
energies,  its  very  good  performance  is  a  confirmation  of  the  reasonableness  of 
its  functional  form.  Two  GGA  functionals  with  an  expansion  in  the  reduced 
gradient  that  is  correct  through  second  order,  namely,  SOGGA  and  SOGGA1 1 , 
were  also  tested  and  found  to  produce  the  best  surface  formation  energies  of  all 
tested  GGA  functionals,  although  they  significantly  overestimate  the  adsorption 
energies. 

We  presented  two  new  exchange-correlation  functionals  for  hybrid  Kohn-Sham 
electronic  structure  calculations  based  on  the  nonseparable  functional  form 


introduced  recently  in  the  N12  and  MN12-L  functionals  but  now  with  the  addition 
of  screened  Hartree-Fock  exchange.  The  first  functional  depends  on  the  density 
and  the  density  gradient  and  is  called  N12-SX;  the  second  functional  depends  on 
the  density,  the  density  gradient,  and  the  kinetic  energy  density  and  is  called 
MN12-SX.  Both  new  functionals  include  a  portion  of  the  Hartree-Fock  exchange 
at  short-range,  but  Hartree-Fock  exchange  is  screened  at  long  range.  The 
accuracies  of  the  two  new  functionals  are  compared  to  those  of  the  recent  N12 
and  MN12-L  local  functionals  to  show  the  effect  of  adding  screened  exchange, 
are  compared  to  the  previously  best  available  screened  exchange  functional, 
HSE06,  and  are  compared  to  the  best  available  global-hybrid  generalized 
gradient  approximation  (GGA)  and  to  a  high-performance  long-range-corrected 
meta-GGA. 

The  new  density  functionals  M11-L  and  M1 1  and  the  SMD  implicit  solvation 
model  were  used  to  predict  reduction  potential-pH  equilibrium  diagrams 
(Pourbaix  diagrams)  for  two  ruthenium-based  water-oxidation  catalysts  and  their 
derivatives  in  aqueous  solution.  The  observed  improved  accuracy  for  transition- 
metal  complexes  opens  new  opportunities  for  the  use  of  theory  in  the 
understanding  and  design  of  catalysts  containing  transition  metals. 

We  reported  a  test  of  30  density  functionals,  including  several  recent  ones,  for 
their  predictions  of  69  singlet-to-singlet  excitation  energies  of  11  molecules.  The 
reference  values  are  experimental  results  collected  by  Caricato  et  al.  for  30 
valence  excitations  and  39  Rydberg  excitations.  All  calculations  employ  time- 
dependent  density  functional  theory  in  the  adiabatic,  linear-response 
approximation.  As  far  as  reasonable,  all  of  the  assignments  are  performed  by 
essentially  the  same  protocol  as  used  by  Caricato  et  al. ,  and  this  allows  us  to 
merge  our  mean  unsigned  errors  (MUEs)  with  the  ones  they  calculated  for  both 
density  functional  and  wave  function  methods.  We  found  21  of  the  30  density 
functionals  calculated  here  have  smaller  MUEs  for  the  30  valence  states  than 
what  they  obtained  (0.47  eV)  for  the  state-of-the-art  EOM-CCSD  wave  function. 

In  contrast,  for  all  of  density  functionals  the  MUE  for  39  Rydberg  states  is  larger 
than  that  (0.1 1  eV)  of  EOM-CCSD.  Merging  the  30  density  functionals  calculated 
here  with  the  26  calculated  by  Caricato  et  al.  makes  a  set  of  56  density 
functionals.  Averaging  the  unsigned  errors  over  both  the  valence  excitations  and 
the  Rydberg  excitations,  none  of  the  56  density  functionals  shows  a  lower  mean 
unsigned  error  than  that  (0.27  eV)  of  EOM-CCSD.  Nevertheless,  two  functionals 
are  successful  in  having  an  overall  mean  unsigned  error  of  0.30  eV,  and  another 
nine  are  moderately  successful  in  having  overall  mean  unsigned  errors  in  the 
range  0.32-0.36  eV.  Successful  or  moderately  successful  density  functionals 
include  seven  hybrid  density  functionals  with  41%  to  54%  Hartree-Fock 
exchange,  and  four  range-separated  hybrid  density  functionals  in  which  the 
percentage  of  Hartree-Fock  exchange  increases  from  0%  to  19%  at  small 
interelectronic  separation  to  65%-100%  at  long  range. 

Using  quantum  chemical  approximations  to  understand  and  predict  complex 
transition  metal  chemistry,  such  as  catalytic  processes  and  materials  properties. 


is  an  important  activity  in  modern  computational  chemistry.  High-level  theory  can 
sometimes  provide  high-precision  benchmarks  for  systems  containing  transition 
metals,  and  these  benchmarks  can  be  used  to  understand  the  reliability  of  less 
expensive  quantum  chemical  approximations  that  are  applicable  to  complex 
systems.  We  studied  the  ionization  potential  energy  of  Fe  and  FeC  and  the  bond 
dissociation  energies  of  FeC  and  FeC"  by  15  density  functional  approximations: 
M05,  M06,  M06-L,  w  B97,  w  B97X,  w  B97X-D,  t  -HCTHhyb,  BLYP,  B3LYP,  M08- 
HX,  M08-SO,  SOGGA1 1 ,  SOGGA1 1-X,  Mil ,  and  Ml  1-L.  All  of  the  functionals 
predict  the  correct  spin  state  as  the  ground  state  of  neutral  iron  atom,  but  five  of 
them  predict  the  wrong  spin  state  for  Fe+  .  In  the  final  analysis,  four  functionals, 
namely  Ml  1-L,  r-HCTHhyb,  SOGGA11,  and  M06-L,  have  small  mean  unsigned 
errors  when  averaged  over  two  bond  dissociation  energies  and  two  ionization 
potentials.  In  fact,  the  results  showed  that  Ml  1-L  gives  the  smallest  averaged 
mean  unsigned  error,  i.e..  Ml  1-L  is  the  most  reliable  density  functional  for  these 
iron  carbide  systems  among  those  studied. 

Kohn-Sham  density  functional  theory  is  in  principle  an  exact  formulation  of 
quantum  mechanical  electronic  structure  theory,  but  in  practice  we  have  to  rely 
on  approximate  exchange-correlation  (xc)  functionals.  The  objective  of  our  work 
has  been  to  design  an  xc  functional  with  broad  accuracy  across  as  wide  an 
expanse  of  chemistry  and  physics  as  possible,  leading — as  a  long-range  goal — 
to  a  functional  with  good  accuracy  for  all  problems,  i.e.  a  universal  functional.  To 
guide  our  path  towards  that  goal  and  to  measure  our  progress,  we  have 
developed — building  on  earlier  work  of  our  group — a  set  of  databases  of 
reference  data  for  a  variety  of  energetic  and  structural  properties  in  chemistry 
and  physics.  These  databases  include  energies  of  molecular  processes,  such  as 
atomization,  complexation,  proton  addition  and  ionization;  they  also  include 
molecular  geometries  and  solid-state  lattice  constants,  chemical  reaction  barrier 
heights,  and  cohesive  energies  and  band  gaps  of  solids.  We  gathered  many  of 
these  databases  into  four  comprehensive  databases,  two  with  384  energetic  data 
for  chemistry  and  solid-state  physics  and  another  two  with  68  structural  data  for 
chemistry  and  solid-state  physics,  and  we  test  two  wave  function  methods  and 
77  density  functionals  (12  Minnesota  meta  functionals  and  65  others)  in  a 
consistent  way  across  this  same  broad  set  of  data.  We  especially  highlighted  the 
Minnesota  density  functionals,  but  the  results  have  broader  implications  in  that 
one  may  see  the  successes  and  failures  of  many  kinds  of  density  functionals 
when  they  are  all  applied  to  the  same  data.  Therefore,  the  results  provide  a 
status  report  on  the  quest  for  a  universal  functional. 

We  presented  a  new  approach  to  calculating  potential  energy  surfaces  for 
photochemical  reactions  by  combining  self-consistent-field  calculations  for  single¬ 
reference  ground  and  excited  states  with  symmetry-corrected  spin-flip 
Tamm-Dancoff  approximation  calculations  for  multireference  electronic  states. 
The  method  is  illustrated  by  an  application  with  the  M05-2X  exchange-correlation 
functional  to  cis-trans  isomerization  of  the  penta-2,4-dieniminium  cation,  which  is 
a  model  (with  three  conjugated  double  bonds)  of  the  protonated  Schiff  base  of 


retinal.  We  found  good  agreement  with  multireference  configuration  interaction- 
plus-quadruples  (MRCISD+Q)  wave  function  calculations  along  three  key  paths 
in  the  strong-interaction  region  of  the  ground  and  first  excited  singlet  states. 

Endohedral  fullerenes  have  great  potential  for  a  variety  of  technological 
applications.  We  considered  B@C6o  and  showed  that  the  amount  of  charge 
transfer  from  the  semimetal  boron  atom  to  the  cage  is  a  strong  function  of  the 
radial  distance  of  the  atom  from  the  center  of  the  fullerene,  and  it  is  controlled  by 
multistate  conical  intersections  whose  associated  ridge  of  avoided  crossings  has 
the  topology  of  a  Euclidean  sphere.  The  potential  energy  surfaces  of  B@C6o  are 
characterized  by  two  kinds  of  local  minima:  those  with  a  boron  atom  located  in 
the  geometric  center  of  the  fullerene,  and  those  with  a  boron  atom  bound  to  the 
fullerene  inner  wall.  At  the  lowest-energy  minimum,  at  the  center,  the  boron  atom 
is  neutral,  whereas  the  transition  to  the  wall  is  accompanied  by  an  electron 
transfer  from  boron  to  the  fullerene  cage.  The  two  kinds  of  minima  are  separated 
by  a  ridge  of  avoided  crossings  that  forms  a  surface  with  a  nearly  spherical 
shape.  The  properties  of  such  systems  may  be  altered  by  controlling  the 
populations  of  the  two  kinds  of  minima,  for  example,  by  application  of  an  external 
field.  Such  switchable  atom-cage  charge  transfer  may  find  applications  in  novel 
molecular  devices. 

The  3d-series  transition  metals  (also  called  the  fourth-period  transition  metals). 

Sc  to  Zn,  are  very  important  in  industry  and  biology,  but  they  provide  unique 
challenges  to  computing  the  electronic  structure  of  their  compounds.  In  order  to 
successfully  describe  the  compounds  by  theory,  one  must  be  able  to  describe 
their  components,  in  particular  the  constituent  atoms  and  cations.  In  order  to 
understand  the  ingredients  required  for  successful  computations  with  density 
functional  theory,  it  is  useful  to  examine  the  performance  of  various  exchange- 
correlation  functionals;  we  did  this  for  4s^  Sd'^  transition-metal  atoms  and  their 
cations.  We  analyzed  the  results  using  three  ways  to  compute  the  energy  of  the 
open-shell  states:  the  direct  variational  method,  the  weighted-averaged  broken 
symmetry  (WABS)  method,  and  a  new  broken-symmetry  method  called  the 
reinterpreted  broken  symmetry  (RBS)  method.  We  found  the  RBS  method  to  be 
comparable  in  accuracy  with  the  WABS  method.  By  examining  the  overall 
accuracy  in  treating  18  multiplicity-changing  excitations  and  10  ionization 
potentials  with  the  RBS  method,  10  functionals  were  found  to  have  a  mean- 
unsigned  error  of  <5  kcal/mol,  with  a)B97XD  topping  the  list.  For  local  density 
functionals,  which  are  more  practical  for  extended  systems,  the  M06-L  functional 
is  the  most  accurate.  And  by  combining  the  results  with  our  previous  studies  of  p- 
block  and  4d-series  elements  as  well  as  databases  for  alkyl  bond  dissociation, 
main-group  atomization  energies,  and  tt-tt  noncovalent  interactions,  we  found 
five  functionals,  namely,  PW6B95,  MPW1B95,  M08-SO,  SOGGA11-X,  and 
MPWB1K,  to  be  highly  recommended.  We  also  studied  the  performance  of 
PW86  and  C09  exchange  functionals,  which  have  drawn  wide  interest  in  recent 
studies  due  to  their  claimed  ability  to  reproduce  Hartree-  Fock  exchange  at  long 
distance.  By  combining  them  with  four  correlation  functionals,  we  found  the 


performance  of  the  resulting  functionals  disappointing  both  for  3d  transition-metal 
chemistry  and  in  broader  tests,  and  thus  we  do  not  recommend  PW86  and  C09 
as  components  of  generalized  gradient  approximations  for  general  application. 

When  the  spins  of  molecular  orbitals  are  allowed  to  be  aligned  with  different 
directions  in  space  rather  than  being  aligned  collinearly,  the  resulting 
noncollinear  spin  orbitals  add  extra  flexibility  to  variational  optimization  of  the 
orbitals,  and  solutions  obtained  with  collinear  spin  orbitals  may  be  unstable  with 
respect  to  becoming  noncollinear  in  the  expanded  variational  space.  We 
explored  whether  and  in  what  way  the  molecular  orbitals  of  the  Kohn-  Sham 
density  functional  theory  become  noncollinear  when  fully  optimized  for  multi¬ 
reference  molecules,  transition  states,  and  reaction  paths.  (I  note  that  a 
noncollinear  determinant  has  intermediate  flexibility  between  a  collinear 
determinant  and  a  linear  combination  of  many  collinear  determinants  with 
completely  independent  coefficients.  However,  the  Kohn-  Sham  method  is 
defined  to  involve  the  variational  optimization  of  a  single  determinant,  and  a 
noncollinear  determinant  represents  the  limit  of  complete  optimization  in  the 
Kohn-  Sham  scheme.)  We  compared  the  results  obtained  with  the  noncollinear 
Kohn-  Sham  (NKS)  scheme  to  those  obtained  with  the  widely  used  unrestricted 
Kohn-  Sham  (UKS)  scheme  for  two  types  of  multi-reference  systems.  For  the 
dissociation  of  the  MnO  and  NiO  transition  metal  oxides,  we  found  UKS  fails  to 
dissociate  to  the  ground  states  of  neutral  atoms,  while  NKS  dissociates  to  the 
correct  limit  and  predicts  potential  energy  curves  that  vary  smoothly  at 
intermediate  bond  lengths.  This  is  due  to  the  instability  of  UKS  solutions  at  large 
bond  distances.  For  barrier  heights  of  03  ,  BeH2  ,  and  H4  ,  NKS  is  shown  to 
stabilize  the  multi-reference  transition  states  by  expanding  the  variational  space. 
Although  the  errors  vary  because  they  are  closely  coupled  with  the  capability  of 
the  employed  exchange-correlation  functionals  in  treating  the  multi- 
configurational  states,  these  findings  demonstrated  that  results  with  collinear  spin 
orbitals  should  be  further  scrutinized,  and  the  future  development  of  exchange- 
correlation  functionals  for  multi-reference  systems  should  incorporate  the 
flexibilities  of  NKS. 

It  is  important  to  test  methods  for  simulating  water,  but  small  water  clusters  for 
which  benchmarks  are  available  are  not  very  representative  of  the  bulk.  We 
presented  benchmark  calculations,  in  particular  CCSD(T)  calculations  at  the 
complete  basis  set  limit,  for  water  26-mers  drawn  from  Monte  Carlo  simulations 
of  bulk  water.  These  clusters  are  large  enough  that  each  water  molecule 
participates  in  2.5  hydrogen  bonds  on  average.  The  electrostatically  embedded 
three-body  approximation  with  CCSD(T)  embedded  dimers  and  trimers 
reproduces  the  relative  binding  energies  of  eight  clusters  with  a  mean  unsigned 
error  (MUE,  kcal  per  mole  of  water  molecules)  of  only  0.009  and  0.015  kcal  for 
relative  and  absolute  binding  energies,  respectively.  Using  only  embedded 
dimers  (electrostatically  embedded  pairwise  approximation)  raises  these  MUEs 
to  0.038  and  0.070  kcal,  and  computing  the  energies  with  the  Mil  exchange- 
correlation  functional,  which  is  very  economical,  yields  errors  of  only  0.029  and 


0.042  kcal. 


In  order  to  understand  what  governs  the  accuracy  of  approximate  exchange- 
correlation  functionals  for  intrinsically  multiconfigurational  systems  containing 
metal  atoms,  the  properties  of  the  ground  electronic  state  of  CaO  werestudied  in 
detail.  We  first  applied  the  T1,  TAE(T),  B1,  and  M  diagnostics  to  CaO  and 
confirmed  that  CaO  is  an  intrinsically  multiconfigurational  system.  Then,  we 
compared  the  bond  dissociation  energies  (BDEs)  of  CaO  as  calculated  by  49 
exchange-  correlation  functionals,  three  exchange-only  functionals,  and  the  HF 
method.  To  analyze  the  error  in  the  BDEs  for  the  various  functionals,  we 
decomposed  each  calculated  BDE  into  four  components,  in  particular  the 
ionization  potential,  the  electron  affinity,  the  atomic  excitation  energy  of  the  metal 
cation  to  prepare  the  valence  state,  and  the  interaction  energy  between  prepared 
states.  We  found  that  the  dominant  error  occurs  in  the  calculated  atomic 
excitation  energy  of  the  cation.  Third,  we  compared  dipole  moments  of  CaO  as 
calculated  by  the  53  methods,  and  we  analyzed  the  dipole  moments  in  terms  of 
partial  atomic  charges  to  understand  the  contribution  of  ionic  bonding  and  how  it 
is  affected  by  errors  in  the  calculated  ionization  potential  of  the  metal  atom.  We 
then  analyzed  the  dipole  moment  in  terms  of  the  charge  distribution  among 
orbitals,  and  we  found  that  the  orbital  charge  distribution  does  not  correlate  well 
with  the  difference  between  the  calculated  ionization  potential  and  electron 
affinity.  Fourth,  we  examined  the  potential  curves  and  internuclear  distance 
dependence  of  the  orbital  energies  of  the  lowest-energy  CaO  singlet  and  triplet 
states  to  analyze  the  near-degeneracy  aspect  of  the  correlation  energy.  The 
most  important  conclusion  is  that  the  error  tends  to  be  dominated  by  the  error  in 
the  relative  energies  of  s  and  d  orbitals  in  Ca+  ,  and  the  most  popular  density 
functionals  predict  this  excitation  energy  poorly.  Thus,  even  if  they  were  to 
predict  the  BDE  reasonably  well,  it  would  be  due  to  cancellation  of  errors.  The 
effect  of  the  cation  excitation  energy  can  be  understood  in  terms  of  an  orbital 
picture,  as  follows.  For  most  functionals  the  predicted  cation  excitation  energy  is 
too  small,  so  it  is  too  easy  to  delocalize  charge  from  the  oxygen  2p  orbital  to  the 
Ca+  d  orbital;  this  overestimates  the  covalency  and  explains  why  most 
functionals  overestimate  the  bond  energy. 

Vanadium  dimer  is  a  notoriously  difficult  case  for  Kohn-  Sham  (KS)  density 
functional  theory  with  currently  available  approximations  to  the  exchange- 
correlation  (xc)  functionals,  and  many  approximate  xc  functionals  yield  an 
exceedingly  large  error  in  the  calculated  bond  energy.  In  this  paper,  we  first  test 
the  bond  energies  estimated  by  43  xc  functionals  and  the  Flartree-  Fock  (FIF) 
method.  The  results  further  confirmed  the  large  errors  and  showed  that,  with  the 
experimental  bond  energy  being  64.2  kcal/mol,  the  KS  calculations  give 
predictions  all  over  the  map  with  errors  ranging  from  -61 .5  to  +60.5  kcal/mol,  and 
the  HF  method  performs  much  worse  with  an  error  of  -124.4  kcal/mol!  The 
reason  for  these  very  large  errors  was  examined  by  analyzing  the  atomic  and 
molecular  orbital  energies  calculated  by  various  xc  functionals.  The  results 
showed  that  the  errors  in  estimates  of  the  bond  energy  of  vanadium  dimer  can 


primarily  be  related  to  the  calculated  energy  gap  between  the  4s  and  3dz2  atomic 
orbitals  of  the  vanadium  atom  and  especially  to  the  3dz2  orbital  energy.  This 
interesting  relation  between  the  errors  in  the  calculated  bond  energy  and  the 
magnitudes  of  the  single-particle  orbital  energies  provides  a  constructive 
alternative  to  the  common  but  more  sterile  explanation  that  it  is  the  static 
correlation  energy  due  to  multicenter  left-right  correlation  that  makes  the 
vanadium  dimer  and  many  other  transition  metal  compounds  so  difficult  for 
Kohn-  Sham  calculations.  One  of  the  most  important  factors  in  determining  the 
critical  atomic  orbital  energy  is  the  amount  of  nonlocal  HF  exchange  that  is 
included  in  the  xc  functional,  but  it  is  still  difficult  to  explain  why  different  local 
functionals  (functionals  with  no  HF  exchange)  yield  quite  different  results.  We 
concluded  that  improving  calculations  of  orbital  energies  of  atoms  may  provide  a 
route  to  improving  the  accuracy  of  theoretical  predictions  of  molecular  bond 
energies  for  systems  containing  metal  atoms. 

We  studied  the  C04O4  subnanocluster  and  its  MeCN-coated  species  using 
density  functional  theory,  and  we  found  that  the  C04O4  core  presents  distinctive 
structures  in  bare  and  ligand-coated  species.  We  proposed  a  possible  ligand- 
mediated  ring  ^  cube  transformation  mechanism  during  the  ligand-coating 
process  of  the  C04O4  core  due  to  the  stronger  binding  energies  of  the  MeCN 
ligands  to  the  3D  distorted  cube  structure  than  to  the  2D  ring  and  ladder 
structures;  theory  indicates  that  three  ligands  are  sufficient  to  stabilize  the  cube 
structure.  Both  ring  and  cube  structures  are  ferromagnetic.  Our  finding  is 
potentially  useful  for  understanding  the  catalysis  mechanism  of  C04O4  species, 
which  have  important  applications  in  solar  energy  conversion  and  water  splitting; 
these  catalysis  reactions  usually  involve  frequent  addition  and  subtraction  of 
various  ligands  and  thus  possibly  involve  core  rearrangement  processes  similar 
to  our  findings. 

We  presented  a  new  theoretical  framework,  called  Multiconfiguration  Pair- 
Density  Functional  Theory  (MCPDFT),  which  combines  multiconfigurational  wave 
functions  with  a  generalization  of  density  functional  theory  (DFT).  A 
multiconfigurational  self-consistent-field  (MCSCF)  wave  function  with  correct  spin 
and  space  symmetry  was  used  to  compute  the  total  electronic  density,  its 
gradient,  the  on-top  pair  density,  and  the  kinetic  and  Coulomb  contributions  to 
the  total  electronic  energy.  We  then  used  a  functional  of  the  total  density,  its 
gradient,  and  the  on-top  pair  density  to  calculate  the  remaining  part  of  the 
energy,  which  we  call  the  on-top  density-functional  energy  in  contrast  to  the 
exchange-correlation  energy  of  Kohn-Sham  DFT.  Because  the  on-top  pair 
density  is  an  element  of  the  two-particle  density  matrix,  this  goes  beyond  the 
Hohenberg-Kohn  theorem  that  refers  only  to  the  one-particle  density.  To 
illustrate  the  theory,  we  obtained  first  approximations  to  the  required  new  type  of 
density  functionals  by  translating  conventional  density  functionals  of  the  spin 
densities  using  a  simple  prescription,  and  we  performed  post-SCF  density 
functional  calculations  using  the  total  density,  density  gradient,  and  on-top  pair 
density  from  the  MCSCF  calculations.  Double  counting  of  dynamic  correlation  or 


exchange  does  not  occur  because  the  MCSCF  energy  is  not  used.  The  theory  is 
illustrated  by  applications  to  the  bond  energies  and  potential  energy  curves  of  H2, 
N2,  F2,  CaO,  Cr2 ,  and  NiCI  and  the  electronic  excitation  energies  of  Be,  C,  N,  N"", 
O,  O'",  Sc"",  Mn,  Co,  Mo,  Ru,  N2,  FICFIO,  C4FI6,  c-CsFle,  and  pyrazine.  The  method 
presented  has  a  computational  cost  and  scaling  similar  to  MCSCF,  but  a 
quantitative  accuracy,  even  with  the  present  first  approximations  to  the  new  types 
of  density  functionals,  that  is  comparable  to  much  more  expensive  multireference 
perturbation  theory  methods. 

Systematic  evaluation  of  the  accuracy  of  exchange-correlation  functionals  is 
essential  to  guide  scientists  in  their  choice  of  an  optimal  method  for  a  given 
problem  when  using  density  functional  theory.  We  determined  the  accuracy  of 
one  Generalized  Gradient  Approximation  (GGA)  functional,  three  meta-GGA 
functionals,  one  Nonseparable  Gradient  Approximation  (NGA)  functional,  one 
meta-NGA,  and  three  hybrid  GGA  functionals  was  evaluated  for  calculations  of 
the  closest  interatomic  distances,  cohesive  energies,  and  bulk  moduli  of  all  3d, 

4d,  and  5d  bulk  transition  metals  that  have  face  centered  cubic  (fee),  hexagonal 
closed  packed  (hep),  or  body  centered  cubic  (bcc)  structures  (a  total  of  27 
cases).  Our  results  show  that  including  the  extra  elements  of  kinetic  energy 
density  and  Hartree-  Fock  exchange  energy  density  into  gradient  approximation 
density  functionals  does  not  usually  improve  them.  Nevertheless,  the  accuracies 
of  the  Tao-Perdew-Staroverov-Scuseria  (TPSS)  and  M06-L  meta-GGAs  and 
the  MN12-L  meta-  NGA  approach  the  accuracy  of  the  Perdew-Burke-Ernzerhof 
(PBE)  GGA,  so  usage  of  these  functionals  may  be  advisable  for  systems 
containing  both  solid-state  transition  metals  and  molecular  species.  The  N12 
NGA  functional  is  also  shown  to  be  almost  as  accurate  as  PBE  for  bulk  transition 
metals,  and  thus  it  could  be  a  good  choice  for  studies  of  catalysis  given  its 
proven  good  performance  for  molecular  species. 


